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Express Mail Label No. EV 396081514 US 
[0001] ADAPTIVE ALGORITHM FOR A CHOLESKY APPROXIMATION 

[0002] CROSS REFERENCE TO RELATED APPLICATION(S) 

[0003] This application is a continuation of U.S. Patent Application Serial No. 

10/153,477, filed May 22, 2002, which is incorporated by reference as if fully set forth. 

[0004] FIELD OF INVENTION 

[0005] This invention generally relates to solving linear systems. In particular, 

the invention relates to using Cholesky approximations to perform a matrix inversion 

of a Hermitian, block-banded and block Toeplitz matrix to solve linear systems. 

[0006] In code division multiple access (CDMA) communication systems, multiple 

communications may be simultaneously sent over a shared frequency spectrum. Each 

communication is distinguished by the code used to transmit the communication. 

[0007] In some CDMA communication systems to better utilize the shared 

spectrum, the spectrum is time divided into frames having a predetermined number of 

time slots, such as fifteen time slots. This type of system is referred to as a hybrid 

CDMA/time division multiple access (TDMA) communication system. One such 

system, which restricts uplink communications and downlink communications to 

particular time slots, is a time division duplex communication (TDD) system. 

[0008] Two approaches to receive the communications transmitted within the 

shared spectrum is single user detection (SUD) and multiuser detection (MUD). SUD 

is used when all received communications experience the same channel response. This 

condition typically occurs in the downlink where all communications are transmitted 
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from the base station or in the uplink where only one user transmits at a particular 

time. In these cases, all the transmitted communications experience a common 

channel response. To compensate for the common channel response, all of the received 

signals are passed through a channel equalization stage. After equalization, the data 

for each communication is recovered using that communications code. 

[0009] The single user data detection problem is typically modeled per Equation 

1 

r = Hs + n 

Equation 1 

r is the received vector. H is the channel response matrix for all the users, n is the 
noise vector, s is the spread data symbols as per Equation 2. 
s = Cd 

Equation 2 

C is the spreading codes of all the users' communications and d is the data vector. 
[0010] Two common approaches to solve Equation 1 are a least squares errors 

(LSE) solution and a minimum mean squares (MMSE) solution. Equation 3 is the LSE 
solution. 



s = (HHH)-iHH £ 

Equation 3 

( ) H is the complex conjugate transpose (Hermitian) operation. Equation 4 is the 
MMSE solution. 
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s = (H H H + a 2 I) 1 H H r 

Equation 4 

o 2 is the noise variance and I is the identity matrix. 

[0011] After solving for s, the data vector d is determined by despreading per 

Equation 5. 

d = O 1 s 

Equation 5 

[0012] In multiuser detection, all the communications' data is determined 

together. Multiuser detection can be used when communications experience differing 
or a common channel response(s). A multiuser detector uses the, known or determined, 
codes of the multiple communications and the determined channel responses to 
estimate the data of all the communications. 

[0013] The multiuser detection problem is typically modeled per Equation 6. 

r = Ad + n 

Equation 6 

A is the system response matrix, which is a convolution of the channel response matrix, 
H, and the code matrix, C. 

[0014] Two common approaches to solve Equation 6 are a zero-forcing (ZF) 

solution and a minimum mean square error MMSE solution. Equation 7 is the ZF 
solution. 

d = (AHA)-iAHr 

Equation 7 
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Equation 8 is the MMSE solution. 

d = (AHA + o2I)-iAHr 

Equation 8 

[0015] A brute force solution to Equations 3 and 4 for SUD or Equations 7 and 8 

have an extremely high complexity. One approach to reduce complexity is Cholesky 
decomposition. Cholesky decomposition is explained in conjunction with the linear 
equation of Equations 9 and 10. 
E = (R)' 1 x 

Equation 9 

x = o H q 

Equation 10 

For a ZF solution, R is o H o. For a MMSE solution, R is o H o + a 2 I. 
[0016] Performing a matrix inversion is a complex operation. To reduce the 

complexity, a Cholesky factor, G, is used when the matrix R is Hermitian block-banded 
and block Toeplitz. G and its conjugate transpose, G 11 , are upper and lower triangular 
matrices. The Cholesky factor, G, is per equation 11. 
GG H = R 

Equation 11 

[0017] After the Cholesky factor, G, is determined, backward substitution is 

performed to determine y per Equation 12. 
y = GHq 

Equation 12 
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[0018] To determine p, forward substitution is performed per Equation 13. 

x = yG 

Equation 13 

[0019] For SUD, p becomes r and q becomes s in Equations 9 and 10. For a ZF 

solution, R is H H H and o H is H H . For a MMSE solution, R is H H H + a 2 I and o H is H H . 
[0020] For MUD, p becomes r and q becomes d in Equations 9 and 10. For a ZF 

solution, R is A H A and o H is A H . For a MMSE solution, o H o is A H A + a 2 I and o H is A H . 
[0021] Although the linear system solution based on Cholesky decomposition 

reduces complexity, determining the Cholesky factor G still involves considerable 
complexity. To reduce the complexity in determining the Cholesky factor G, an 
approximate Cholesky factor G is used. Approaches to determine G have been to 
determine a column or row of G and to replicate that column or row to generate the G 
matrix. The approximation of G may introduce an error into the data detection 
process. This error may reduce performance of the receiver to an unacceptable level. 
[0022] Accordingly, it is desirable to have alternate approaches for determining 

the Cholesky factor. 

[0023] BACKGROUND 
[0024] 

[0025] SUMMARY 

[0026] Data is to be detected in a wireless communication system. A plurality of 
communication signals are received. A solution for estimating data of the received 
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communication signals is modeled using a linear system requiring a matrix inversion. 
Columns or rows of an approximate Cholesky factor are determined. A difference 
between the determined columns or rows is determined. If the determined difference is 
less than a threshold, subsequent columns or rows are determined by previously 
determined columns or rows. The data of the received communication signals is 
estimated using the approximate Cholesky factor. 
[0027] 

[0028] BRIEF DESCRIPTION OF THE DRAWING(S) 

[0029] Figure 1 is an embodiment of a communication system using an 

approximate Cholesky factor. 

[0030] Figure 2 is a flow chart of determining an approximate lower triangular 

Cholesky factor by columns. 

[0031] Figure 3 is an illustration of an approximate lower triangular Cholesky 

factor determined by columns. 

[0032] Figure 4 is an illustration of an approximate lower triangular Cholesky 

factor determined by rows. 

[0033] Figure 5 is an illustration of an approximate upper triangular Cholesky 

factor determined by columns. 

[0034] Figure 6 is an illustration of an approximate upper triangular Cholesky 

factor determined by rows. 

[0035] Figure 7 is a flow chart of generally determining a Cholesky factor either 
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by columns or rows. 

[0036] DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 
[0037] The preferred embodiments of the invention can generally be used with 

any type of CDMA system, such as a TDD/CDMA, TDMA/CDMA or frequency division 

duplex/CDMA communication system, as well as other types of communication 

systems. The determination of the Cholesky factor can be used to perform a matrix 

inversion in a linear system for a variety of applications. 

[0038] Figure 1 illustrates an embodiment of a communication system using an 

approximate Cholesky factor. A transmitter 20 and a receiver 22 are shown in Figure 
1. The transmitter 20 may be located at a user equipment or multiple transmitting 
circuits 20 may be located at the base station. The receiver 22 may be located at either 
the user equipment, base station or both. 

[0039] Data symbols to be transmitted to the receiver 22 are processed by a 

modulation and spreading device 24 at the transmitter 20. The spreading and 
modulation device 24 spreads the data with the code(s) assigned to the 
communication(s) carrying the data. The communication(s) are radiated by an 
antenna 26 or antenna array of the transmitter 20 through a wireless radio interface 
28. 

[0040] At the receiver 22, the communication(s), possibly along with other 

transmitters' communications, are received at an antenna 30 or antenna array of the 
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receiver 22. The received signal is sampled by a sampling device 32, such as at the 
chip rate or at a multiple of the chip rate, to produce the received vector, r. The 
received vector r is processed by a channel estimation device 36 to estimate the channel 
impulse responses, H, for the received communications. The channel estimation device 
36 may use a training sequence in the received communication, a pilot signal or 
another technique to estimate the impulse responses. A data detection device using an 
approximate Cholesky factor 34, such as a SUD or MUD, uses the code(s) of the 
received communication(s) and the estimated impulse response(s) to estimate soft 
symbols, d, of the spread data. 

[0041] For a data detector using SUD, the data detector uses the received vector 

r and channel response matrix H (and for a MMSE solution the noise variance o 2 ) to 
determine the spread data, s_using an approximate Cholesky factor. The spread data 
vector s is despread using the user codes C. For a data detector 34 using MUD, the 
data detector 34 uses the received vector r and the system response matrix A (and for a 
MMSE solution the noise variance a 2 ) to determine the spread data, dousing an 
approximate Cholesky factor. The system response matrix A is determined using the 
channel responses H and the user codes C. 

[0042] The data detector 34 uses an approximate Cholesky factor to estimate the 

data as per Figure 2. The approximate Cholesky factor is either an upper or a lower 

triangular matrix. The following discussion addresses one embodiment which repeats 

block columns in a lower triangular G factor to produce an approximate Cholesky 

factor. However, the same approach can be used to repeat rows in a lower triangular 
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matrix or repeat columns or rows in an upper triangular matrix. 
[0043] Figure 3 illustrates an approximate Cholesky factor being lower 

triangular and produced by repeating columns. The Cholesky factor matrix has blocks 
of dimension K by K. K is the number of user communications being processed. Each 
block column of the matrix has a predetermined length, such as L or 2L. L is the 
number of symbols that influence each other by intersymbol interference (ISI). 
Preferably, each block column is 2L blocks in length, which simulations indicate has a 
negligible degradation depending on the fading channel. 

[0044] The factor is preferable calculated starting with the leftmost column, step 

38. Each block of the column is determined, preferably, from top block to bottom block. 
After the first column is determined, the determination of blocks in the next column is 
started, preferably from top to bottom. As shown, the next column is one block lower in 
the matrix. Entries of the leftmost block column are compared to entries of the next 
column, such as per a threshold test, step 40. If the comparison of the two columns 
indicates that the differences between the columns is less than the threshold, the 
remaining columns are produced by repeating and moving down one block the last 
produced column, step 44. If the comparison is greater than the threshold, the next 
column is calculated, step 46 and 40. This just calculated column is compared to the 
previous column. This process is repeated until the two last produced columns pass the 
comparison test or an upper limit in the produced columns is reached. If the upper 
limit is reached, the remaining columns are produced by repeating the last fully 
produced or upper limit column, step 48. 
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[0045] Although the entire columns may be compared or a sampling of entries in 

each column, preferably, the top blocks in each column are compared. By comparing 
only the top blocks, the rest of the column can be produced by replication of the 
previous column. 

[0046] One threshold test calculates an error between the top block of the two 

columns, column A and B. The error is determined by an error function, which takes 

the difference between the normalized A and B block matrices. One example of an 

error function is per Equationl4. 

/ A norm(A - B) 

error(A, B) = . - 

^norm(A) • norm(B) 

Equation 14 

The notation norm denotes the matrix norm. If the error is less than the threshold, the 
remaining columns are produced by repeating the last fully determined column. If the 
error exceeds the threshold, the rest of the column is determined along with the top 
block of the next column for use in another threshold test. 

[0047] Using the threshold test allows the Cholesky factor to be determined with 

a certain accuracy. As a result, the accuracy of the data detection can similarly be 
controlled. 

[0048] Although computing the Cholesky factor was illustrated for repeating 

columns for a lower triangular matrix, the repeating can be performed for a row or 
column of a lower or upper triangular matrix. Figure 4 illustrates repeating the rows 
of a lower triangular matrix, Figure 5 illustrates repeating columns of an upper 

-10- 



1-2-0276. 1US 



triangular matrix and Figure 6 illustrates repeating the rows of an upper triangular 
matrix. Since the first rows of a lower triangular matrix and the first columns of an 
upper triangular matrix do not have all the elements of a representative row/column, 
the threshold test is not performed until the L th row/column is determined. 
[0049] In general, the determining of the Cholesky factor is as follows as shown 

in Figure 7 for either a row or column of a lower or upper triangular matrix. Blocks of 
the matrix are determined from a corner along the main diagonal. Although it is 
preferred to use the upper left most corner, the lower right most corner may be used. 
Blocks of the first row/column, at the corner, are determined, step 50. Sequentially, 
the blocks of the neighboring rows/columns are determined, step 52. Once a 
row/column of full length is determined, a comparison test is performed on each 
subsequently produced row/column, step 54. If the two rows/columns are similar 
enough, subsequent rows/columns are produced by repeating the previous row/column, 
step 56. If the comparison is not favorable, a next row/column is determined, step 58. 
This procedure is repeated until a favorable comparison results or, alternately, an 
upper limit is reached. If the an upper limit is reached, subsequent rows/columns are 
produced by repeating the last completed row/column, step 60. 

The approximate Cholesky factor is used to determine the data of the received 
communications, such as in a SUD or MUD using ZF or MMSE data detection 
approaches. 

* * * 



